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O ' Abstract 

^■"■^ ' We present results concerning the non-perturbative evaluation of the ghost-gluon running 

QCD coupling constant from Nj = 2 twisted-mass lattice calculations. A novel method for 
calibrating the lattice spacing, independent of the string tension and hadron spectrum is 

r> I presented with results in agreement with previous estimates. The value of A.j^ is computed 

j^ ' from the running of the QCD coupling only after extrapolating to zero dynamical quark mass 

and after removing a non-perturbative OPE contribution that is assumed to be dominated 
by the dimension- two (A^) gluon condensate. The effect due to the dynamical quark mass 
in the determination of Aj^ is discussed. 

LPT-Orsay 10-37 
UHU-FT/10-31 



^Unite Mixte de Recherche 8627 du Centre National dc la Recherche Scientifique 



Contents 

1 Introduction 

2 The running coupling in Taylor scheme 

2.1 Pure perturbation theory 

2.2 OPE power corrections 

3 The lattice computation of the Taylor coupling 

3.1 The lattice action 

3.2 The computation of the gluon and ghost Green functions . 

3.3 On the treatement of lattice artefacts 

3.3.1 Hypercubic i?(4)-extrapolation 

3.3.2 Quark mass artefacts 

4 Computing Aj^ and the gluon condensate 

4.1 Looking for the ^^ plateau" 

4.2 Global fit and the calibration of lattice spacing 

4.3 The contribution from the Wilson coefficient higher orders 

4.4 Discussing the systematical uncertainties 

4.4.1 Volume effects 

4.4.2 Three-loop versus four-loop confrontation 

4.4.3 The impact of higher-orders in the OPE expansion 

4.5 Conversion to physical units and quark mass effect . . . . 



5 Discussions and conclusions 25 



A Appendix: The Wilson coefficients at the four-loops order |28 



1 Introduction 

QCD is believed to be the theory of the strong interactions with, as only inputs, one mass pa- 
rameter for each quark species and the value of the QCD coupling constant at some energy or 
momentum scale in some renormalization scheme (Alternatively, this last free parameter of the 
theory can be fixed by Aqcd, the energy scale used as the typical boundary condition for the 
integration of the Renormalization Group equation for the strong coupling constant). This is the 
parameter which expresses the scale of strong interactions, the only parameter in the limit of 
massless quarks. While the evolution of the coupling with the momentum scale is determined by 
the quantum corrections induced by the renormalization of the bare coupling and can be computed 
in perturbation theory, the strength itself of the interaction, given at any scale by the value of 



the renormalized coupling at this scale, or equivalently by Aqcd, is one of the above mentioned 
parameters of the theory and has to be taken from experiment. 

The QCD running coupling can be also obtained from lattice computations, the free parameters 
being adjusted from experimental numbers, masses, decay constants etc. These parameters being 
settled, the lattice calculation of Aqcd proceeds in several manners: the implemention of the 
Schrodiger functional scheme (see, for instance, [I]-|l] and references therein), those based on 
the perturbative analysis of short-distance sensitive lattice observables as the "boosted" lattice 
coupling (see for instance [5H8] and reference therein) and, in particular, those based on the study 
of the momentum behaviour of Green functions(see P4T6] and references therein) are among 
the most extendedly applied. Indeed, the confrontation of the behaviour with respect to the 
renormalization scale of 2-gluon and 3-gluon Green functions with the corresponding perturbative 
predictions leaves us with a good estimate of as-, its running leading to the determination of 
Aqcd, but also reveals a dimension- two non-zero gluon condensate in the Landau gauge. The 
possible phenomenological implications in the gauge-invariant world of such a dimension-two gluon 
condensate and in connection with confinement scenarios has been also largely investigated (see 
for instance [T7]). 

In [18], the Green's function approach proposed in ref. [16] was followed exploiting a non- 
perturbative definition of the coupling derived from the ghost and gluon propagators for the 
determination of hrg^ in pure Yang-Mills [Nf = 0). In that work, the renormalization scheme for 
the ghost-gluon vertex corresponding to the latter coupling was properly defined. The quenched 
lattice results were analyzed over a wide momentum window, applying a ^^plateau" procedure to 
extract simultaneously both Ay[^ and the gluon condensate. The result is consistent with other 
calculations, for the description of the gluon and ghost Green functions and for the running of the 
strong coupling in ^^quenched^ QCD. 

In the present paper we extend [18] to the case in which twisted Nf = 2 dynamical quarks are 
included in the lattice simulations for several different bare lattice couplings {(3 = 3.9,4.05,4.2) 
and different dynamical quark masses. We use the configurations produced by the ETM Col- 
laboration [19] . This offers the opportunity to study the effect of the quark mass on the lattice 
determination of the strong coupling, and we will see that this effect if far from negligible. Similar 
works have started some time ago [15] using unquenched lattice configurations with, at first, rather 
heavy Nf = 2 dynamical quarks, and continuing more recently [31[IB] in a more realistic case. The 
comparison of our current results with those needs certainly due account of the dependence on 
the dynamical masses. 

Let us now summarize our strategy. For every value fi of the momentum scale we compute, 
from the lattice simulations, the value of the strong coupling constant. This can be converted via 
a four loops formula to a value for Aj^(yU,). Ay^ is a scale independent constant which sets the 
strong interaction scale. It results that Aj^ (/i) should be independent of fi as soon as we are in the 
perturbative regime. As we shall see, this is far from true at energies of several GeV's, which are 
generally believed to lie in the perturbative regime. This surprising feature was already noticed in 
the quenched case [18]. We then need to take into account non-perturbative contributions using 



Wilson expansion. In Landau gauge there exists only one dimension two operator : A^ = A'^A'Ij^. 
The Wilson coefficient of that operator has been computed to order a^ [20]. We will assume 
that only this A^ operator contributes q We fit < A^ > so that, once the non-perturbative 
contribution subtracted, one gets a good "plateau" for A^jg(yu). We thus get an estimate of both 
A^jg and < A^ >. 

The paper is organized as follows. In section [21 we outline and discuss all the analytical 
tools, perturbation theory and Wilson expansion, needed to describe the running coupling in the 
appropriate renormalization scheme. In section [3l we give the details of the lattice computation 
of the coupling, describe the treatment of lattice artefacts and depict the analysis procedure 
leading to the estimate of Aj^ and the gluon condensate. The analysis is performed in section 
m where we also present the results and discuss different sources of systematical uncertainties 
(special attention is paid to the higher order contribution in both the perturbative and the OPE 
expansions) . Finally, we conclude in section |5l 

2 The running coupling in Taylor scheme 

Among the many possibilities to the compute a strong coupling as from lattice simulations, it 
has been shown [18] that the so-called Taylor scheme is among the most tractable ones because, 
with the help of the so-called Non-renormalization Taylor theorem, the coupling can be computed 
from two-point Green functions renormalized in MOM scheme. Following the usual notation we 
will write Landau gauge gluon and ghost propagators as: 

(G'^'):(/.A) - ^*4*..-^ 

(fWj-^/.A) = -i.,:^^; (2.i; 



with A the regularisation cutoff (a ^(/3) if, for instance, we specialise to lattice regularisation). 
The renormalized dressing functions, Gr and Fr are defined through : 

Gn{p',fi') = limZ3-l(/i^A)G(p^A) 

A— >CXD 

Fnip^^i') = limZ3-l(/x^A)F(p^A), (2.2) 

with standard MOM renormalization condition 

G'fi(/i2,/i2) = F^(/i^/i2) = l. (2.3) 

Then, we can consider the ghost-gluon vertex which could be non-perturbatively obtained through 
a three-point Green function, defined by two ghost and one gluon fields, with amputated legs after 



It is indeed not easy with present accuracy to discard higher order operators. 



dividing by two ghost and one gluon propagators. This vertex can be written quite generally as: 

a q-k 
rf (-g, k-q-k)= -^Q^- = ig.r^'^ {q.H,{q, k) + {q - k),H2{q, k)) , (2.4) 

k q 

where Hi^2 are the two involved independent scalar form factors, go is the bare strong coupling, q 
is the outgoing ghost momentum and k the incoming one. This bare three-pointt Green function 
will be renormalized according to: 

Tr = Z^V. (2.5) 

In the MOM renormalization procedure, as explained in ref. [IB], Zi is fully determined by de- 
manding that one specific combination of those two form factors (chosen at one's will) be equal 
to its tree-level value for a specific kinematical configuration; while 

,«(/.^) = hm Z-^(/.^A^),o(A^) = hm ^3^'(/^^ A^)^3y , A^) ^^^^.^ ^^.6) 

A^oo 3 A^o^ Zi(/i2,A2) 

where one puts explicitly the cut-off dependence. If one turns now to the Taylor scheme, i.e. 
a specific MOM-type renormalization scheme defined by a kinematical configuration with zero 
incoming ghost momentum, one obtains [21] 



Zi(/i2,A2) =1; (2.7) 

and one is left with 

«t(^^) = ^ = hm ^G(/.^ k')F\^^\ A^) ; (2.8) 

where we also apply the renormalization condition for the propagators, eqs. (I2.2|2.3p . to replace 
the renormalization constants, Z^ and Z^, by the bare dressing functions. As emphasized in [IS] , 
the remarkable feature of eq. (12. 8p is that it involves only F and G so that no measure of the 
ghost-gluon vertex is needed for the determination of the coupling constant. 

2.1 Pure perturbation theory 

The Taylor coupling and the one renormalized in the standard MS prescription, as any other 
different definitions of the coupling constant, can be related through: 

«T(/i^) = «(/i^) fl + E^^(^)^) ' (2-9) 



where the coefficients q's can be obtained in perturbation theory [22|l23] 
507 - 40iV/ 
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It was proven in ref. [IB] that these coefficients could be also directly derived from the anomalous 
dimensions for gluon and ghost propagators, as eq. (12. 8p indicates. The three coefficients in 
eq. (12.1 op obviously define unambigously the running of ar up to four-loops given by 
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with t = \ia.jT, since the coefficients of the /3- function of a^, 



i=0 



can be derived, 



/3o = /3o = ll-^iV^ 
A = ^i = 102-yiV; 

= 3040.48 - 625.387 Nf + 19.3833 A^J 

/33 = ^3-2^2Cl+^lC?+^o(2C3-6C2Ci+4c?) 

= 100541 - 24423.3 Nf + 1625.4 Nf - 27.493 A^J , (2.13) 

from the knowledge of those coefficients, q's, and from that of the standard MS /3- function, 

^ i=o ^ ^ 
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given at four loops in ref. [21] {(3q and /3i being scheme-independent). As for the Aqcd parameters 
in the two schemes, they are related through 

ci 507 - AONf 

^ = e"2^ = e~792-48iV/ ^ 0.541449 . (2.15) 

Then, the lattice data for the coupling, obtained through eq. (12. Sp . can be confronted to the 
perturbative formulae, eq. (12. lip with the /3- function coefficients given by eq. (I2.10p . over the 
large-momentum window where the four-loop perturbation theory is reliable. 

2.2 OPE power corrections 

As was previously done for the quenched analysis in [18], in order to take full advantage of the 
lattice data (and reduce the systematic uncertainties) when confronting them with a formula 
for the QCD running coupling, one needs to take into account the gauge-dependent OPE power 
corrections (cf. [T2HI1]) to ot- In Landau gauge there exists only one dimension- two operator 
allowed to have a vacuum expectation value: A^ = A^ = Ai^AI^^. We will stick to it, leaving the 
discussion of higher dimension operators to section HJ We will also not go beyond the leading logs 
in this section, leaving again the discussion of higher orders [20] to section HI 

The leading power contribution to the ghost and gluon propagators can thus be computed 
using the operator product expansion [25] (OPE), as is done in ref. ^E\, and one obtains: 

{G^'Xu^q') = {Gf^XM') + < A{^^) + ■ ■ ■ ^2.16) 

where the Wilson coefficients, diagramatically expressed as follows (the bubble means contracting 
the color and Lorentz indices of the incoming legs with l/25st5ar) 



w""^ 





Ku = ^QOOO^fiji^MIOOQ^ + 2x ^qq!?OOQCLQC&)Q^ ' (217) 



can be computed by invoking the SVZ factorisation [27] . Thus, after some algebra and the 
appropriate renormalization at the subtraction point q^ = fi"^, according to the MOM scheme 



definition (details of tlie computation can be found in [l3l[T8l|26]), one gets at tree level 

G„(,^^,^) = G^,„i„^p.') (i + i|M|hkij , (2J8) 

where the multiplicative correction to the purely perturbative contributions were determined up 
to corrections of the order l/g'^ or lng//i. Finally, putting together the defining relation eq. fl2.8p 
and the results eqs. (12.1 8p we obtain 

arifi') = lim ^^^^(/.^ A)G(/x^ A) 

A-i-oo 47r 

= hm ^F2(go^A)F^(/i^go2)G(go^A)G«(/i^go') 



„pcrt^^2^ i;,2 / 2 2^ ^ / 2 2x f i , ^ 9Ti%yA ) R,ql 






where ^q ^ Aq^^q is some perturbative scale and the /3- function, and its coefficients in eq. (I2.13p . 
of course describe the running of the perturbative part of the evolution, a^^^. The anomalous 
dimension for the Wilson coefficient, 

^A^iai^^')) = hm -^lnZ^.(/i^ A^) = -^f ^ +••• (2-20) 

A^oo dm/i^ Air 

where Aj^ = Z~^J^, is neglected in eq. (12.191) . 

The leading logarithm contribution for the Wilson coefficient are incorporated as explained 
in [18], yielding: 




oF-\qi)) 4(iV2 



c 



where 7^ can be taken from [2ni[2H] to give 



7^ 105-8iVf 9 
1 - ^ = 1 = 5 , 2.22 

which agrees for Nf = with the power of the logarithmic correction applied, and shown to have 
a negligible impact on a in [18]. We shall first apply formula eq. (I2.2ip . approximated up to the 
four-loop level in perturbation and up to the leading-log in the OPE expansion, to describe the 
lattice data in the next sections. As already mentioned we postpone the use of four-loop Wilson 
coefficients [20] and the study of the impact of higher order operators to the end of section HI 
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3 The lattice computation of the Taylor coupling 

The following section is devoted to the computation of the running coupling in Taylor scheme, 
eq. f|2.8p . from the lattice. The results presented here are based on the gauge field configurations 
generated by the European Twisted Mass Collaboration (ETMC) with the tree-level improved 
Symanzik gauge action [29] and the twisted mass fermionic action [30] at maximal twist. 

3.1 The lattice action 

A very detailed discussion about the twisted mass and tree-level improved Symanzik gauge actions, 
and about the way they are implemented by ETMC, can be found in refs. [T9 |[3TH33] . Here, for the 
sake of completeness, we will present a brief reminder of the twisted action and the run parameters 
for the gauge configurations that will be exploited in the present work (See tab. [T]). 

The Wilson twisted mass fermionic lattice action for two flavours of mass degenerate quarks, 
reads (in the so called twisted basis 



(3.1) 



X 

D^ = -7, (V, + v;) - y v,v; , 

where ttiq is the bare untwisted quark mass and /Xg the bare twisted quark mass, t^ is the third 
Pauli matrix acting in flavour space and r is the Wilson parameter, which is set to r = 1 in the 
simulations. The operators V^ and V* stand for the gauge covariant nearest neighbour forward 
and backward lattice derivatives. The bare quark mass rriQ is related as usual to the so-called 
hopping parameter /t, by k = 1/(8 + 2amo)- Twisted mass fermions are said to be at maximal 
twist if the bare untwisted mass is tuned to its critical value, merit- This is in practice done by 
setting the so-called untwisted PCAC mass to zero. 

In the gauge sector the tree-level Symanzik improved gauge action (tlSym) [29] is applied. 
This action includes besides the plaquette term U^/^^j^ also rectangular (1 x 2) Wilson loops U^^^^,. 
It reads 

^^ = f E ( ^0 E il - »eTr(f/2,^^:j}+6i X; {1 - MeTr(f/i;iy } j , (3.2) 

X \ ^,U=1 IJ,,U = 1 / 

where /3 = Q/qq, go being the bare lattice coupling and it is set bi = —1/12 (with bo = 1 — 8bi as 
dictated by the requirement of continuum limit normalization). Note that at 6i = this action 
becomes the usual Wilson plaquette gauge action. The run parameters for (3 and /ig of the gauge 
configurations that will be exploited in the following can be found in tab. [T] 



/3 


afiq 


Volume 


Number of confs. 


3.9 


0.004 

0.0064 

0.010 


243 X 48 


120 

20 

20 


4.05 


0.003 
0.006 
0.008 
0.012 


32^ X 64 


20 
20 
20 
20 


4.2 


0.0065 


323 X 64 


200 



Table 1: Run parameters of the exploited data from ETMC collaboration. 
3.2 The computation of the gluon and ghost Green functions 



Applying eq. f l2.8p demands to compute the gauge-fixed 2-point gluon and ghost Green functions 
from the lattice. To this goal, we exploited ETMC gauge configurations obtained for f3 = 3.9, 
P = 4.05 and /3 = 4.2 and a large variety of dynamical quark masses, fixed by the values of the Hq 
parameter. The lattice gauge configurations are transformed to Landau gauge by minimising the 
following functional of the SU(3) matrices, U^{x), 



Fu[g] = Re 



E E Tr (1 - ^9i^)U,{x)g\x + /i)^ 



(3.3) 



with respect to the gauge transform g, by applying a combination of overrelaxation algorithm and 
Fourier acceleration Q This procedure does not avoid the possibility of lattice Gribov copies that, 
in any case, have been reported to have a nonsignificant influence beyond the lowest momenta. 
Then, the gauge field is defined as 



A^{x + fi/2) 



U,{x) 



Uji^) _ 1 ^^ f U,{x)-Ul{x) 
2iago 3 I 'iiago 



(3.4) 



where fi indicates the unit lattice vector in the direction /i and go is the bare coupling constant. 
The 2-gluon Green functions is computed in momentum space by 



where 



) indicates the Monte-Carlo average and where 
1. 



Kip) 



-Tr 



J2 M^ + A/2) exp{ip{x + ii/2))X' 



(3.5) 



(3.6) 



•^We end when |(9^jj4^p < 10 ^^ and when the spatial integral of Aq is constant in time to better than 10 ^ . 
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A'* being the Gell-Mann matrices and the trace being taken in the 3x3 color space. 
On the other hand, the ghost propagator is also computed in Landau gauge, 

{F^'r\^-y) - i{M-X)^ (3-7) 

as the inverse of the Faddeev-Popov operator, that is written as the lattice divergence, 

M(f/) = -Iv ■ 5(t/) (3.8) 

where the operator D acting on an arbitrary element of the Lie algebra, t] reads: 

D{U)v{x) = ^ {U^{x)r]{x + /i) - v{x)U^{x) + r/(x + fi)Ul - Ul{xUx)) . (3.9) 

More details on the lattice procedure for the inversion of Faddeev-Popov operator can be found 
in 135). 



3.3 On the treatement of lattice artefacts 

The lattice estimates of the Landau-gauge propagators through eqs. (I3.5f3.7p . after Fourier trans- 
forming the ghost correlator and appropriate projection of both as indicated by eq. (12. ip . lead 
to the determination of ghost and gluon dressing functions to be used in eq. (12. 8p . Both dress- 
ing functions are dimensionless lattice correlation functions (let us note both as Q = F, G) that, 
because of general dimensional arguments, depend on the lattice momentum ap^, where 

27Tn 
Pm = ^ n = 0,l,---,N, (3.10) 

and on the strong interaction scale Aqcd- Anticipating on the averaging over hypercubic orbits 
and on the treatment of hypercubic lattice artefacts we get Q = Q{a? p^, a^Aq^j-,). Our choice for 
the lattice action ensures that the discretization artefacts due to the lattice are 0{a?), where a is 
the lattice spacing. 

Thus, the running coupling in Taylor scheme is obtained by 

aT{^^') = lml^F2(a^p^a^A^eD)G(a^P^a^A^CD) , (3.11) 

where taking the limit of a vanishing lattice spacing indeed implies the proper elimination of 
lattice artefacts. 
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3.3.1 Hypercubic iJ(4)-extrapolation 

A first kind of artefacts tliat can be systematically cured [361137] are those due to the breaking 
of the rotational symmetry of the Euclidean space-time when using an hypercubic lattice, where 
this symmetry is restricted to the discrete H{4) isometry group. It is convenient to compute first 
the average of any dimensionless lattice quantity Q{ap^) over every orbit of the group H{4). In 
general several orbits of H{A) correspond to one value of p^. Defining the ^^(4) invariants 

pw=x:p; p'''=i:pi (3-12) 

it happens that the orbits of H(4:) are labelled u by the set p^, a^p^^', a^p^^'. In the continuum limit 
the effect of a^pW, a^p'^^ vanishes. We can thus define the quantity Q{ap^) averaged over -^(4) as 

g(a^ p^aV^aV^a^ A^cd)- (3.13) 

If the lattice spacing is small enough such that e = c?'p^ j-j? <^ 1, the dimensionless lattice 
correlation function defined in eq. fl3.13p can be expanded in powers of e: 

Q{a'p\ aV'Ka%^'\ a^A^co) = Q(«V, ^^A^cd) + ^ a'^ + ■ ■ ■ (3.14) 

ae ^^g p 

if (4) methods are based on the appearance of a O(a^) corrections driven by ap'^^ term. The basic 
method is to fit from the whole set of orbits sharing the same p^ the coefficient dQ/de and get the 
extrapolated value of Q, free from -ff (4) artefacts. If we further assume that the coefficient 

K{a p , a AqcdJ = 



de 

has a smooth dependence on a^p"^ over a given momentum window, we can expand R as R = 
Rq + RiO?p^ and make a global fit in a momentum window between (p — 5, p + 5) to extract the 
extrapolated value of Q for the momenta p in the window, and shift to the next window etc. This 
procedure of fitting with windows is somehow different from the basic one, since the extrapolation 
does not rely on any particular assumption for the functional form of R. On the other, the 
systematic error coming from the extrapolation can be estimated by modifying the width of the 
fitting window. 

It is worthwile to mention that we considered in this work anisotropic lattice of the type L^xT, 
with T = 2L. This finite volume effect reduces the H{4) lattice symmetry to H{3). Deviations 
from H[4:) are to be expected in the long-distance physics. But ultraviolet physics should not be 
affected. As far as we are interested in the high-momentum regime, we will assume the previous 
treatement of the lattice artefacts to be valid. 



^On totally general grounds, any iJ(4)-invariant polynome can be written only in terms of the four invariants 
p[^'l with i — 1, 2, 3, 4 [3nil3Z]- As a consequence of the upper cut for momenta, the first three of these invariants 
suffice to label all the orbits we deal with and hence any presumed dependence on p^^' is neglected. 
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3.3.2 Quark mass artefacts 



In this section we will consider the influence of the dynamical quark masses. We will argue that 
this is a 0{a^fi^) effect and therefore that it is a lattice artefact. 

Following the treatement described in the previous section, we have calculated the if (4)-free 
ghost and gluon dressing functions, that we shall denote in the following by F, G. These dressing 
functions can be combined in order to calculate the if(4)-free lattice coupling through eq. (13. lip . 

In the analysis performed ref. [18] by exploiting "quenched" configurations, a pretty good 
scaling of cty was found, computed at different /3's from eq. (13. lip , once one parameter describing 
the lattice spacing ratios had been fitted. This seems to imply that the residual O{o?'p'^) artefacts 
are negligible for the "quenched" coupling constant in Taylor scheme after if (4)-extrapolation has 
been performed. Of course 0{o?h^Qj^) artefacts may still be hidden in the matching coefficient 
between lattice spacings. 

It is worth pointing that all the divergent contributions appearing for the dressing functions 
or the bare coupling when the lattice spacing vanishes cancel when combined in eq. (13.111) . 



|3=3.9 



P=4.05 



H 0.45 




a(p)p 




Figure 1: The Taylor couplings estimates through eq. ()3.1ip . after iJ(4)-extrapolation, at /3 = 3.9 for 
^lq = 0.004, 0.0064, 0.010 (left-hand plots) and at /3 = 4.05 for /ig = 0.003, 0.006, 0.008, 0.012 (right-hand 
plots). 

However, when analyzing "unquenched" lattice configurations, one should keep in mind that 
one additional mass scale, the dynamical quark mass, is playing a role. In Fig. [1] one can see the 
Taylor coupling after hypercubic extrapolation for different /ig at fixed /9 = 3.9 and 4.05. Indeed, 
a dependence in /ig is clearly seen. If it is an artefact the dependence should be in a?^i^q. If it is 
an effect in the continuum it should be some unknown function of the physical mass /Zg. Trying 



13 



an 0(a yUq) dependence, we write the expansion : 



^ / 2 2 2 2\ 

aria p ,a fi ) 



9M 



rii 2 2 2 2\ 172/ 2 2 2 2\ 

G{a p ,a /ig)F [a p ,a /ij 



aT(aV,0) + 



dai 



dia^fi'] 



/ 2 2\ 2 2 I 

[a p ) a i^g + 



(3.15) 



tr,,»L 



a{P)p= I. OS 




(afP) U )- 



a(|3)p=l.36 














0..M 








- 


• a(|i)p = 


us 


0J35 






0.315 


^ 


"~T~~^^ 


i- 


_: 



X 10" 0.^65 



(a(P) K,)" 





(ii(|3) U ) 



a<|!)p=1.52 




wpj^r 



(a(P) U )" 



(a(P) U )" 



Figure 2: We plot the values of the Taylor coupling at /? = 4.05, computed for some representative values 
of the lattice momentum, a(4.05)p = 1.08, 1.18, 1.24, 1.36, 1.45, 1.52, in terms of a^(4.05)/Xg and show the 
suggested linear extrapolation at a^/i^ = 0. 



Provided that the first-order expansion in eq. (I3.15P is reliable, a linear behaviour on a^/i^ has 
to be expected for the lattice estimates of Sy for any fixed lattice momentum computed from 
simulations at any given /3 and several values of /ig. We explicitely check this linear behaviour to 
occur for the results from our P = 4.05 and /3 = 3.9 simulations and show in Fig. [2]some plots of St 
computed at /3 = 4.05 (where four different quark masses are available) for some representatives 
lattice momenta in terms of a^yU^. We thus write the Taylor expansion as and after neglecting the 
0{a'^) contributions get 



^/22 2 2\ /2\iD/2 2\ 22 

aria p ,a ft ) = ar(p ) + Ro{a P ) a n 



where Ro{a?p^) is defined as 



Ro{aV) ^ 



dciT 



(3.16) 



(3.17) 
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In fig. [31 we plot Ro{a'^p'^) as a function of ap computed for the four lattices simulations at /3 = 4.05 
with different quark masses and for the three ones at /3 = 3.9 (see tab. [1]). Indeed, it can be seen 
that a constant behaviour appears to be achieved for p > pmin — 2.8 GeV. We will not risk an 
interpretation of the data below (ap)min- The striking observation here is that above Pmin both 
lattice spacings exhibit a fairly constant Rq{o?p'^) and a good enough scaling between both /3's. 



500 






-500 




0.5 



a(f3)p 



Figure 3: The slopes for the mass squared extrapolation in terms of ap computed for the four lattices 
simulations at /3 = 4.05 (32^ x 64) with a/ig = 0.003, 0.006, 0.008, 0.012 and for the three ones at /3 = 3.9 
(24^ X 48) with a/ig = 0.004, 0.0064, 0.010. The supposed constant behaviour appear to be reached when 
lattice volume effects become negligible for each simulation. 

The analysis of the slopes (see fig. [3]) leaves us with both a fair estimate of Rq above 2.8 
GeV, Rq ~ —0.9 ■ 10^, and an intrinsic definition for the momentum window where eq. (I3.16P can 
be applied to extrapolate at any momentum. Thus, after the extrapolation with the previously 
obtained Rq down to vanishing a^q, one obtains the three estimates for the running coupling at 
/3 = 3.9, 4.05, 4.2 plotted in terms of the momentum in lattice units, ap, in fig. HI (a), which shows a 
very smooth running behaviour. These are the lattice estimates for the coupling to be confronted 
to the analytical prediction given by eq. 02.2ip . 

The fact that Rq with our present data goes to the same constant for both /3's, leads us to 
consider that the /i, dependence of a is mainly a lattice artefact (else it should be a function of 
/ig and not of a^q). The slope Rq is not small. We did not expect this. It has to be seriously 
considered as it affects the result on Aj^. The slope being negative, the extrapolation to vanishing 
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afiq leads to a larger value for Aj^ than if we had estimated it at finite a fig. We shall see that 



this effect is of the order of an increase of 40 MeV on A- 
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Figure 4: Left-hand plot (up): The Taylor coupling, free of H{4) and mass-quarks artefacts, computed 
by applying eq. (j3.16p for the three (3 = 3.9, 4.05, 4.2 and plotted in terms of the lattice momentum a{/3)p. 
The three other plots show Atto in lattice units, computed by the inversion of eq. 



with the lattice 



couplings plotted in the upper left-hand plot, for /3 = 4.2 (upper right), f3 = 4.05 (bottom left) and 
/3 = 3.9 (bottom right); the black circles are for a perturbative inversion (with zero gluon condensate) 
and blue squares are computed with the best-fit of the gluon condensate (see the text). 

Taking into account the effects due to dynamical quarks in a global analysis of the lattice 
determinations is among the main results of this paper. This will lead to a proper extrapolation to 
the continuum limit. 
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4 Computing Aj^ and the gluon condensate 

In this section, following [18], we will apply a "^/atean" -procedure exploiting eq. f l2.2ip to get a 
reliable estimate of the AqcD-pararneter from the Taylor running coupling constant computed as 
explained in the previous section from lattice simulations with Nf = 2 dynamical quark flavors. 

4.1 Looking for the ^'plateau'^ 

In fig. m we show the estimates of A^^ obtained when interpreting the lattice coupling computed 
from eq. f l3.1ip for /3 = 3.9,4.05,4.2 for all momenta inside the window where the slope for the 
/ig-extrapolation behaves as a constant, i.e. for p > pmin ~ 2.8 GeV, up to our chosen lattice upper 
bound Ifl of ap < 1.6. The estimate of Aj^ is done first (black circles) thanks to the inverted four- 
loop perturbative formula for the coupling, eq. (12. lip . These estimates systematically decrease as 
the lattice momentum increases, while if we were in the perturbative region it should be a constant 
as Ay^ is, by definition, a constant in the perturbative expansion. This clearly reveals the necessity 
of applying the non-perturbative formula including power corrections, eq. (I2.2ip . with a non-zero 
gluon condensate. This is also done in fig. HI where the same is plotted but inverting instead the 
non-perturbative formula (blue squares). The value of the gluon condensate has been determined 
by requiring a "p/ateaw" to exist over the total momenta window. More precisely, one requires the 
best-fit to a constant of the estimates of Aj^, in lattice units, in terms of the lattice momentum, 

(xi^Vi) = iapi,Ai) , (4.1) 

where, of course, i runs to cover all the lattice estimates of the coupling inside the defined window 
and where Aj is obtained by inverting U 

«r6og^l = ^ii^!4^, (4.2) 



A? 



a2p2 



with a"^"^ given by the perturbative four-loop formula eq. (12. lip . Q'5"(a^p^,0) taken from the 
extrapolation of the lattice couplings by eq. (I3.16P and c resulting from the best-fit (it appeared 
written in terms of the gluon condensate in eq. (I2.19P ) of Aj to a constant. Thus, Aj is required 
to reach a "p/ateau" , behaving as a constant when i runs, in terms of the lattice momentum. 

This procedure leaves us with estimates of Ay^ and the gluon condensate (computed from 
the best-fit determination of c), expressed in lattice units, for any /3. However, we will take this 
as a striking illustration of the necessity of including non-perturbative power corrections in the 
analysis (see fig. H]) but we will only report the results obtained in the next section when a global 
fitting strategy will be applied to our /Xg-extrapolated lattice data for the three different /3's. 



^Above some value of ap the lattice artefacts become overwhelming. One must choose an upper bound, this 
choice being to some extent arbitrary. 

^For the sake of simplification we are using here the tree level value for the Wilson coefficient of the dimension-2 
condensate, i.e. a constant ; higher orders will be considered in the next sections. 
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4.2 Global fit and the calibration of lattice spacing 

The running of ar given by the combination of Green functions in eq. (13. lip and the extrapolation 
through eq. fl3.16p . provided that we are not far from the continuum hmit and discretization errors 
are treated properly, depend only on the momentum (except, maybe, finite volume errors at low 
momenta). The supposed scaling of the Taylor coupling implies for the three curves plotted in 
fig- 111(a) to match to each other after the appropriate conversion of the momentum (in x-axis) 
from lattice to physical units, with the multiplication by the lattice spacing at each /3. Thus, we 
can apply the "plateau" -method described in the previous subsection for the three /3's all at once 
by requiring the minimisation of the total x^: 

f a(B) a(3)\ ' U^ - ^AM^mX 
«(/i^o)AMS.c,;;7^,;;7^] - 2^ 2^ ^^77^"^ , (4-3) 

j=0 



X >— -Ms,^,^(^^),^(^^)y ^Z^ ^2(A^ 



where the sum over j covers the sets of coupling estimates for the three /3's (/3o = 3.9, /3i = 4.05, 
P2 = 4.2), the index i runs to cover the fitting window defined, as previously explained, through 
the slope analysis |j and Aj(/3j) is again obtained for any (3j by requiring 



port / 1 '^ {Pj)Pi 

"t log 



A?(/3.) 



/ / /^ xN 2 /io„-^!l^M\ "A 



1 + 



<^mv r^'^^^mr 



^ a-iPMUMj l^log^^^ J 



aT{a\(3,)plO) (4.4) 



where now we apply the OPE formula including the leading logarithm correction for the Wilson 
coefficient, eq. (I2.2ip . with a^'^ given by the perturbative four-loop formula, eq. (12. lip , and where 
0'{(3o)qo = 4.5 (this means go ~ 10 GeV). The errors for the extrapolated couplings are estimated 
by jackknife analysis and properly propagated through the perturbative inversion to give S{Ai). 
The function x^ is minimised over the functional space defined by the four parameters that are 
explicitly put in arguments for eq. (I4.3p 's l.h.s.: a(/3o)Aj^, c, ^|j4, ^W)- Thus we obtain all at 
once Aj^ and the gluon condensate, in units of the lattice spacing for /3o = 3.9, and the ratios of 
lattice spacings for our three simulations after the extrapolation to the limit /ig — )■ (see tab. [2]). 
The errors are calculated again by jackknife analysis. 

The ratios of lattice spacings can be applied to express the momenta for all the three sets 
of coupling estimates plotted in fig. H] (upper left-handed plot) in units of the lattice spacing at 
f3 = 3.9. Thus they indeed match each other and fit pretty well to the analytical prediction given 
by eq. (I2.2ip with the best-fit parameters for Aj^ and the gluon condensate, in units of l/a(3.9) 



'^In the case of /32 = 4.2, as only the simulation for one quark mass, /i^ = 0.0065, is exploited, one extrapolates 
by applying the slope Rq computed for /3o = 3.9 and /3i = 4.05 all the coupling estimates obtained inside the same 
lattice-momentum window determined for /3i = 4.05. So we do because simulations at both /3i = 4.05 and /32 — 4.2 
where performed in 32'^ x 64 lattices and the impact of volume effects were supposed to determine the lower bound 
of the fitting window. 
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Figure 5: The scaling of tlie Taylor coupling computed by applying eq. ()3.16p for the three /3 = 
3.9, 4.05, 4.2 is shown. The lattice momentum, a(/3)p in the x-axis, is converted to a physical momentum 
in units (the same for the three /3's) of a(3.9)~^. Fo so to do, the ratios of the lattice spacings must be ap- 
plied. The solid curve is the non-perturbative prediction given by eq. (j2.2ip with the best-fit parameters 
for Aj^ and the gluon condensate, and the dotted one is the same but with zero gluon condensate. 



(see tab. [2]), as can be seen in the plot of fig. [5] . The quality for the fits drastically deteriorates as 
data below a(3.9)p ~ 1.2 are included, whichever value results for the gluon condensate. Thus, a 
fitting window excluding those data are applied in obtaining the best-fit parameters in tab. [2] and 
the best-fit curve in fig. [5] which clearly detaches from the data below such a lower bound. 



4.3 The contribution from the Wilson coefficient higher orders 

The Wilson coefficients for gluon and ghost propagators have been very recently obtained at four- 
loop level [20] . As it is shown in appendix, by exploiting the results of this ref . [2U] , the four-loop 
OPE formula for the T-scheme coupling for Nf = 2 can be obtained by replacing in eq. (12.211) : 



9t{(1o) 



-)■ 



X 



aUQo) (l + 1-2932 a(/i2) + 1.9363 a^(/i^) + 3.8296 a^(/i^)) 
(1 - 0.7033 a(g2) _ 0.3652 a'^{q^) + 0.0051 a^{ql)) 



(4.5) 
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This paper 


String tension 


a(3.9)/a(4.05) 


1.224(23) 


1.255(42) 


a(3.9)/a(4.2) 


1.510(32) 


1.558(52) 


a(4.05)/a(4.2) 


1.233(25) 


1.241(39) 


Ams«(3-9) 


0.134(7) 




^2(A2)a2(3.9) 


0.70(23) 





Table 2: Best-fit parameters for the ratios of lattice spacings, Aj^ and the gluon condensate (for which 
a(3.9)(7o = 4.5 is chosen). For the sake of comparison, we also quote the results from [38] that were 
obtained by computing the hadronic quantity, ro/a(/3), and applying to it a chiral extrapolation. 

where a = a"^"^ and Qq is the renormalization momentum for the local operator A^ (see next 



eq. ( ]A.9|) and compare with eq. (I2.2ip that only incorporates the leading logarithmic contribution) 
which we fixed, as previously indicated, by requiring: a(3.9)go = 4.5 . 

Then, we can repeat the analysis of previous sections, after the replacement prescribed in 
eq. (14. 5p . and obtain the results collected in tab. |3l Thus, a strong stability results from the 
fits for the estimates of Aj^ (all of them being compatible within the statistical uncertainties and 
varying less than a 2.3 %) and a fairly convergent behaviour for that of the gluon condensate which, 
computed with the one-loop Wilson coefficient, clearly borrows something from next-to-leading 
contributions. 





One loop 


Two loops 


Three loops 


Four loops 


Ams«(3.9) 


0.134(7) 


0.136(7) 


0.137(7) 


0.138(7) 


g'{A')a'{3.9) 


0.70(23) 


0.52(18) 


0.44(14) 


0.39(14) 



Table 3: Best-fit parameters for Ay^ and the gluon condensate (for which a(3.9)(7o = 4.5 is chosen) 
by applying a OPE formula including the logarithmic corrections for the Wilson coefficient at one, two, 
three and four loop order. No difference is seen for the ratios of lattice spacings. 



4.4 Discussing the systematical uncertainties 

The main sources of systematical uncertainties affecting the determination of the best-fit pa- 
rameters (Aj;jg, the non-perturbative gluon condensate and the ratios of lattice spacings) of the 
matching previously described are expected to come from the truncation of the perturbative series 
for the theoretical prediction of the coupling in eq. (I2.2ip . the possible effect of higher-orders in 
the OPE expansion and from the finite volume effects in lattice simulations. We will pay attention 
in the following to these error sources. 
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4.4.1 Volume effects 

As can be seen in tab. [H we exploited lattice simulations in volumes 24^ x 48 at /3 = 3.9 and 32^ x 64 
at /3 = 4.05 and 4.2. In the case of /3 = 4.2, in order to spare computing time, we use a volume 
smaller than what is usually needed to measure hadronic quantities (in particular, a 48"^ x 96 lattice 
at /3 = 4.2 and afig = 0.002 is required), relying on the hope that, being interested in ultraviolet 
quantities, the finite volume effects will be reduced. One expects that the smaller the product of 
momentum and lattice size, the larger is the volume-effect impact. Indeed, we introduced a lattice 
momentum cut, a(/3)pmin; when we studied the quark- mass extrapolation, which the slopes for 
the mass squared extrapolation detached below from the constant behaviour, and we interpreted 
this as a possible volume effect. Furthermore, in fig. no impact of any remaining volume effect 
on the determination of the coupling is seen: the impressive scaling shown by the results from 
our three simulations in fig |5] seems to confirm that we are finally left with no important volume 
effect. 

4.4.2 Three-loop versus four-loop confrontation 

A standard way to estimate the effect of perturbative-series truncation is to repeat the analysis 
described in previous sections but applying instead a three-loop formula for the perturbative 
inversion. If this is done, one obtains the results collected in tab. HI 





Four loops 


Three loops 


a(3.9)/a(4.05) 


1.224(23) 


1.229(23) 


a(3.9)/a(4.2) 


1.510(32) 


1.510(29) 


a(4.05)/a(4.2) 


1.233(26) 


1.234(25) 


Ams«(3-9) 


0.134(7) 


0.125(6) 


g'{A')a\3.9) 


0.70(23) 


0.80(20) 



Table 4: Best-fit parameters for the ratios of lattice spacings, Aj^ and the gluon condensate (for which 
a{3.9)qo = 4.5 is chosen). 

Then, we can conclude that no noticeable impact from the perturbative truncation is resulting 
on the determination of the ratios of lattice spacing. Nevertheless, a systematic uncertainty of 
roughly a 7 % can be estimated from the discrepancy of the three and four-loops estimates for the 
Aj^. Analogously, the determination of the gluon condensate is affected by a correction of the 
order of 13 %. 

4.4.3 The impact of higher-orders in the OPE expansion 

We previously paid attention to the comparison of the best-fit results when applying Wilson 
coefficients for the OPE expansion at different loop-orders (see tab. [3]). We then concluded that 
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the different estimates of Aj^ differ from each other less than a 2.3 %, while those for the gluon 
condensate fairly converge as the loop-order increases to roughly one half of the one-loop result. 
On the other hand, as previously explained, the quality for the fits drastically deteriorates 
as data below a(3.9)j9 ^1.2 are considered (the x^ becomes of the order of four times larger). 
This could be an indication of the impact of OPE higher power corrections that could be simply 
parameterized as 

ttr,P4(Ai^) = aT(/i^) + ^ , (4.6) 

where ar is given by eq. (I2.2ip and C4 is a constant encoding all the information coming from the 
condensates of higher- dimension operators and the Wilson coefficients (their anomalous dimension 
is thus neglected). If we try a fit with eq. (14. 6 p to the lattice data, good-quality fits are obtained 
for negative values of C4, of the order of —0.1, while the positive estimated contribution of the 
gluon condensate increases drastically. The fitting function with or without the C4///'' turn out to 
be very close to one another over the whole momentum window. In other words there is a valley 
in the parameter space in which both the 1///^ and the l//i^ coefficients vary in an anticorrelated 
way without any significant change of the value of fitting function. Furthermore, including the 
04//^^ term increases drastically the errors. Aj^ is also correlated with the other two parameters: 
it decreases when 04//!^ decreases below zero, by about 10 % from the fit without C4/yU.^ term and 
the fit with it. 

Our conclusion is that the fit with 04//!^ is extremely unstable. We therefore decide not to 
include it in our fits. The resulting systematic uncertainties is not larger than ~ 10% on Aj^. 
The estimate of the gluon condensate may be more severely affected. But, as shown in the 
quenched case [18], the estimates of the gluon condensate stemming from different quantities, 
when neglecting C4//i^ terms, are quite compatible. This would not be possible if the dimension- 
four operators, with different Wilson coefficients for every quantity, were playing a significant role. 
We remain thus rather confident in our estimate oi < A^ >. What should be stressed here is that 
the necessity of a positive non perturbative contribution is unavoidable, and that it is well taken 
into account by the dominant dimension-two < A^ > condensate. 

4.5 Conversion to physical units and quark mass effect 

In this section, we will apply the results for the lattice spacings at /3 = 3.9, 4.05 and 4.2 obtained 
in ref. [3H] through a very exhaustive investigation of the light meson physics using maximally 
twisted mass fermions for Nf = 2 degenerated quark flavours (the gauge configurations we use in 
this work were part of the data ensembles generated by ETMC and analysed in ref. [38]). The 
physical scale is given by requiring /^ = 130.7 MeV as also done in [31]. We first computed the 
ratios of those lattice spacings obtained in ref. [5B] and compared them with the ones obtained in 
this paper in tab. |2J The agreement is indeed remarkable. Then, by applying the result 



a(3.9) = 0.0801(14) fm , (4.7) 
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to convert into physical units our estimates of Ay^ and the gluon condensate (see tab. [2]), the final 
results are: 

/oo^ , ^o , ^^ N /0.0801 fm- 130.7 MeV \ ,, ,, 
A^ = (330 ± 23 ± 22_33) x ( -^^^^^ j MeV , 

9, 9, ,,9, / +7N /0.0801 fm- 130.7 MeV \^ ^ , , , 

/(go)(^\o = (4.2±1.5±0.7+^)x (^ ^^3 g^ ^^ j GeV^ (4.8) 

where we first quote the purely statistical errors, where the one from the lattice size determination 
in ref. [5B] has been properly into account, and then the main systematical uncertainties which 
were detailed in the previous subsection and that, as explained, mainly tend to increase the 
gluon condensate and to decrease the value of Aj^. In particular, for the gluon condensate, 
the contribution to the uncertainty of higher powers in OPE expansion, although unequivocally 
increasing its size, is very hard to be estimated (as discussed in the previous section) and we 
explicitely indicated this by the addition of a question mark in the upper systematical errors of 
eq. (14. 8p . It should be also remembered that go = 4.5 a(3.9)~^ = 11.1 GeV. In ref. [38], the 
lattice spacings result from combined fits including several lattice simulations at different /3's, 
after chiral extrapolations on the quark mass, where the physical scale is fixed by /^ = 130.7 
MeV. In particular, eq. (14. 7p and the ratios we presented for the sake of comparison in tab. [21 
were obtained in that ref. |3B] through a combined fit including /3 = 3.9,4.05 and 4.2, but results 
obtained through combined fits including either /3 = 3.8, 3.9 or /3 = 3.8, 3.9 and 4.05 were also 
reported. Had we applied instead of eq. (14. 7p the other results reported in [38] in order to convert 
our best-fit parameters into physical units, we would obtain the results of tab. [51 Thus, we can 
roughly estimate how our final result in eq. (14. 8 p is systematically affected by the conversion to 
physical units and consider the central value of A^^ roughly to range from 313 to 335 MeV and 
that of the gluon condensate from 3.8 to 4.3 GeV^. Furthermore, in order to make easier any 
further comparison (avoiding also the ambiguities related to the conversion to physical units), 
we collect in tab. [61 the ratios of A^jg and some momentum-dimension physical quantities also 
obtained in ref. [38] • In particular, if we compare our Nf = 2 estimate for Aj^, converted to 
physical units, in eq. (14. 8 p and the same obtained by applying the Schodinger functional method 
in ref. [3], Aj^ = 245(16) (16) MeV, they clearly differ. Nevertheless, had we compared their 
estimate of ^^oA^jg = 0.62(4) with ours in tab. [6l ^-oAy^ = 0.72(5), we would conclude that they 
almost agree with each other within their statistical error intervals. Thus, this indicates that the 
main source of discrepancy for these two results comes from setting the physical scale. 

On the other hand, for the purpose of illustrating the effects derived from the dynamical quark 
mass on the determination of A^jg, we can analyze instead of the /Xg-extrapolated data for the 
coupling (plotted in fig. ^ the ones obtained from the the lattices at /3 = 4.2 with afiq = 0.0065, 
P = 4.05 with afiq = 0.008 and /3 = 3.9 with afiq = 0.010. In view of the results of the lattice 
spacings ratios in tab. [2 and eq. (14. 7p . these assumed to be independent of the quark mass, and 
after the appropriate renormalization, 

f^R = y^ , (4.9) 

Zp{qo) 
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/3's in fits 


a(3.9) (fm) 


Ams (MeV) 


g^{A^) (GeV^) 


3.9, 4.05, 4.2 


0.0801(14) 


330(23) 


4.2(1.5) 


3.9, 4.05 
3.8, 3.9, 4.05 


0.0790(26) 
0.0847(15) 


335(28) 
313(22) 


4.3(1.7) 
3.8(1.4) 



Table 5: Best-fit parameters for Aj^ and the gluon condensate (for wliich a(3.9)go = 4.5 is cliosen) 
by applying the different lattice spacings for a(3.9) obtained in ref. [38] after combined fits including 
simulations with the different /3's indicated in the first column. 





h 


/o 


1/ro 


T^u^d 


A-MS 


2.52(18) 


2.71(19) 


0.72(5) 


92(9) 



Table 6: dimensionless ratios of At 



ijyjg and some momentum dimension quantities taken from ref. [38], 
where the gauge configurations analyzed in this paper were also exploited. Each number in the table is 
obtained by dividing the quantity indicated in the row label (Aj^) over the one indicated in the column 
label. 



where we apply the MS renormalization constant, Zp, at the renormalization scale Qq = 2 GeV 
given in ref. [3S], three very similar quark masses will be obtained for the three simulations (see 
tab. [7]). Then, by implementing the "plateau method" with the OPE formula in eq. (I2.2ip . but 
applying eq. (14. 7p to convert into physical units, one will be left with an estimate of Arg^ at a 
renormalized quark mass of the order of /i/j ~ 50 MeV (see fig. |6]). The result of this analysis is: 



A- 



MS 



294 ± lOMeV 



(4.10) 



that appears to indicate a trend that should be kept in mind to compare with previous (or future) 
unqueched estimates of Aj^: the larger is the quark mass, the lower is the estimate of Aj^. For 



instance, in ref. [15] , Aj^ = 264 MeV was reported as the result of a preliminary analysis of the 
three-gluon coupling from lattice simulations with two fiavours of Wilson dynamical quarks and 
with a renormalized sea-quark mass roughly ranging from 100 to 400 MeV. Both this preliminary 
results and that of eq. (I4.10p seem to be in the right ballpark. So more if one considers the 
uncertainty derived from the physical lattice calibration previously discussed (see tab. E]) and 
which makes the central value of Aj^ to range from 278 to 298 MeV. 

This trend for the behaviour of the estimate for Aj^ from the lattice, in particular from the 
lattice strong coupling in Taylor scheme, is corroborated by the results obtained from the analysis 
of quenched lattice simulations in ref. [TS] (where Aj^ = 224 MeV). This can be clearly seen in the 
plots of fig. El where the OPE formulae with the values of Aj^ and the gluon condensate obtained 
by fits of results from lattices at ^xr = 0, /ij^j ~ 50 in this paper and by a fit of results from quenched 
lattices in ref. [18] appear displayed. It should be noted that the quenching approximation can 
be understood as equivalent to consider infinite dynamical quark masses. It may be also worth to 
remember that, in the three cases, the OPE formulae fit pretty well to the lattice results for the 
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/3 


afXg 


fiR (MeV) 


4.2 


0.0065 


53.4 


4.05 


0.008 


54.9 


3.9 


0.010 


52.8 



Table 7: Quark masses for the three simulations which we re-analyze the results from in order to get an 
indication of the quark mass effect. 

Taylor coupling from momenta roughly ranging from 3 to 6 GeV. 



tt^ 0.5 r 



' 


• n=0 




-- g"<A->=5.9 GeV, 


A,j5=294 McV 


— g^<A">=4.5 GeV, 


A^j=330 MeV 


■■■■ g^<A'>=5.1 GeV, 


A^j=224 MeV 




tt.p 0.5 r 



Hu - 50 McV 
— g"<A'>=4.5 GeV, Ayj=330MeV 
-- g"<A'>=5.9 GeV, Aj,j=294MeV 
■ ■ ■ ■ g"<A'>=5. 1 GeV, A =224 MeV 




^i (GeV) 



|i (GeV) 



Figure 6: The Taylor coupling after the chiral extrapolation plotted in fig. [5l but now in terms of the 
momentum in physical units (left), and the same Taylor coupling computed for the three lattice data 
sets of tab. [3 without chiral extrapolation, also in terms of the momentum in physical units (right). The 
OPE formula with the best-fit parameters obtained for extrapolated lattice data at /x/j = (solid red 
line), for the lattice data at /i/j ^^ 50 MeV (black dashed line) and for the quenched {fig — > oo) lattice 
data in ref. [TB] (blue dotted line) are displayed in both plots for the sake of comparison. 



5 Discussions and conclusions 

We analyzed gauge configurations with twisted Nf = 2 dynamical quarks at the maximal twist 
produced by the ETM collaboration for lattices at several /3's (3.9, 4.05, 4.2), and with several 
different bare twisted masses, and computed the strong coupling constant renormalized in the 
so-called Taylor scheme. The main advantage of this scheme being not to require any three-point 
Green function computation, statistical fluctuations became under enough control to permit a very 
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elaborate treatment of the lattice artefacts and a precise estimate of the couplings at the infinite 
cut-off limit. The coupling estimates for lattices at different /3's were seen to match pretty well, as 
should happen if the cut-off limit is properly taken, when plotted in terms of the renormalization 
momenta converted to the same units by applying the appropriate lattice spacings ratios. These 
ratios could be either taken from independent computations or obtained by requiring the best 
matching with pretty compatible results. Indeed, to require the best matching for the strong 
coupling computed in terms of the momentum for lattices at different /3's was proposed as an 
alternative method to determine the lattice spacings ratios in ref. [18] and shown here to work 
pretty well. 

Thus, once we are left with the estimates of the coupling constant extrapolated at vanishing 
dynamical mass fig, for every value of the renormalization momentum, /i, they were converted 
via a fit with a four loops formula into the value of Aj^. As also noticed for the analysis of 
lattice results with Nf = 0, the values of Aj^ so obtained for every value of fi resulted not to 
be independent on /i over a momentum window roughly ranging from 2 to 6 GeV. This implies 
a non-negligible impact of non-perturbative contributions which we accounted for by a Wilson 
OPE expansion assumed to be dominated by the condensate of the dimension two operator: 
A^. Thus, after including the gluon condensate contribution we recovered a "plateau" for the 
fitted values of Ajj^ in terms of the renormalization momentum and estimated the condensate. 
After discussing the several main sources of systematic uncertainties, specially the higher order 
contributions in both perturbative and OPE expansions, we converted our fitted parameters into 
physical units by applying [3B]: a(3.9) = 0.0801(14) fm. Thus we obtain 330(23) (22) .33 MeVs for 
Aj^ and 4.2(1.5) (0.7) GeV^ for the gluon condensate. As was proven in the quenched case [18], a 
positive dimension two gluon condensate is unequivocally needed to account for the momentum 
behaviour of the strong coupling computed from the lattice. Whether condensates of higher order 
operators have to be included or not in the Wilson expansion is a different matter that cannot be 
properly addressed with our current data. Indeed, a negative contribution of the order 1/p'^ can 
borrow something from the 1/p^ gluon condensate, both OPE formulae with and without such a 
contribution remaining almost totally indistinguishable within our fit momentum window. Then, 
since the fits including 1/p^ contributions become unstable and very noisy, we do not take them 
into account. A possible way-out could be to perform a cross-check by analyzing the running 
behaviour of the renormalization constant of other operators. In particular for the quenched case, 
apart from gluon and ghost operators, the vector part of the quark propagator is studied in ref. [3S] 
and the fitted gluon condensate is found to be compatible with the one obtained from the ghost- 
gluon vertex [18]. A program to study systematically the non-perturbative OPE contributions to 
the renormalization constant of quark operators is in progress [IQ] . 

Furthermore, we also paid attention to the effect of the dynamical quark masses. As we per- 
formed an extrapolation to zero bare twisted mass for the lattice estimates of the strong coupling 
constant, the A^^ result corresponds to a world with Nf = 2 chiral quarks. Nevertheless, for the 
sake of comparison, we applied the same analysis procedure to three sets of lattice configurations 
corresponding to practically the same renormalized dynamical quark mass (roughly 50 MeVs), 
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with no mass extrapolation, and estimated a lighter value of Aj^ = 294(10) MeV (the central 
value can also range from 278 to 298 MeV, depending on which value we use for the lattice spac- 
ing in physical units). This is compatible with previous preliminary results |j [15l[T6]. We may 
recall here that the quenched results [18] for A^jg are still lower, ranging mainly around 230 MeV. 
Since quenching can be understood as the situation with infinite dynamical masses, one may infer 
a general trend that A^^ increases when the quark mass decreases, with a finite value both at 
infinite and vanishing quark mass. 

As a matter of fact, most of the results for Aj^ obtained with Wilson fermions (around 260-270 
MeV; see for instance |5] and references therein) lie below our zero- mass result, but also below 
the phenomenogical value which could be obtained, after the appropriated conversion, from the 
experimental world average, a{Mz) = 0.01184(7) [H]. The standard procedure to convert A^^ to 
Oy^ at a given scale, typically the mass of Z boson, implies the RGE four-loop evolution of the 
coupling and the three-loop matching at the quark thresholds. Provided that the MS running mass 
of the charm quark is 1.5 GeV, the conversion from A^^" to A^~ implies to evolve a^^ over 

an energy range where perturbation theory could fail; but the conversion from Aj-=i.~ to A^^" 
is however out of the scope of the above-described standard conversion procedure. This is why, 
at present , we prefer not to compare with phenomenological results. This task is to be properly 
acomplished when Nf = 2 -|- 1 + 1 and Nf = 4 lattice simulations will be available. Neverteless, 
our current results for Aj^ with two light quarks flavours, after chiral extrapolation, seems to 
point that systematic effects due to the dynamical quark masses could explain the discrepancy 
between the lattice estimates for the strong coupling with Wilson fermions and its experimental 
determination. 
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^The authors of ref. [16] also exploit strong coupling computed on the lattice and renormalized in Taylor scheme, 
although they presented this scheme as a combination of MOM prescription for the propagators and the MS for 
the ghost-gluon vertex. 
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A Appendix: The Wilson coefficients at the four-loops 
order 

The purpose of this appendix is to exploit the four-loops results in ref. [2D] to derive the Wilson 
coefficients with the appropriate renormalization prescription and modify properly eq. fl2.2ip . 
Following [131 HH] the equations fl2.16p for ghost and gluon propagators can be rewritten after 
renormalization as 



GR{q,^l) = Co(-^,a(/i)J+C2(^-^,a(/i,^ 4(iV2 - l)g2 



/i2' '^7 '\fl^' '"^'J 4(iV2_l)^2 



With the help of the appropriate renormalization constants, one can also write eq. (lA.ip in terms 
of bare quantities: 



G(g^A^) = Z3(/i^A^)co(^,a(/i 



^^2'"VA- 7 4(iV2_l)g2 ' 



+ Zsifi',A')Z-Xf,',A') c^l^^.aifi')) ,,Jf\,^, , (A.2) 



where Aj^ = Zj^A?. A totally analogous equation for the ghost dressing function Fi^q^^h?), with 
Cj and Z3 in place of Cj and Z3. Now, as the /^-dependence of both l.h.s. and r.h.s. of eq. (IA.2p 
should match each other for any g, one can take the logarithmic derivative with respect to /x and 
infinite cut-off limit, term by term, on r.h.s. and obtains: 



7(«(/i'))+ TTT^ + /3(«(/^'))^ lncoP^,«(/i^)) = 



-7A<«(/i2))+7K/i^)) + <|TTT^ + /3(«(/^'))|-l Inc2f4.«(/^')) = 0, (A.3) 



dlogji^ da) \/i2 

where 7(a(/i2)) is the gluon propagator anomalous dimension and 7A2(tt(/^^)) is the anomalous 
dimension for the local operator A^ defined in eq. fl2.20p and that was obtained at four-loop in 
ref. [28]. Both eqs. ( lA.Sp can be finally combined to give: 

-7A^(«(/^')) + ^r^ + /5(«(/^'))^ ^, f = , (A.4) 

^log/x2 da] co(j;,a(/i2)) 

and we can proceed in the same way for the ghost dressing function and derive analogous equations 
for the Wilson coefficients, Cj, that differ from those for Cj only because 7(a(/i2)) ^g^j^gg ^^j^g place of 
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7(«(/i^))- Thus, the combination C2/C0 obeys exactly the same eq. (]A.4p . above derived for C2/C0, 
that can be solved by applying the following ansatz, 



C2 



/^^ 



r,a(/^^) 



Co , 2 



,a{iJ^) 



C2 (l,a(g^)) 






i 
\ i / 



(A.5) 



where we use that the leading Wilson coefficient is to be renormalized in the MOM renormalization 
prescription such that Co(l, a;(/i^)) = 1 and where the exponent a and the coefficients rj's are 
required to satisfy eq. (IA.4p . Concerning the boundary condition for C2, the prescription applied 
for the renormalization of the local operator A^ to obtain its four-loops anomalous dimension, 
7^2 , in [28] is the standard MS and, with this prescription, C2 is computed at the four- loop order 
in ref. [20] (see the eq. (7) of that paper). We only need to take g^ = /i^ in the expression given 
in ref. [20] and have thus 02(1, a(/i^)). Then, one obtains at the four-loop order: 



C2 l^^,a(/i2; 
Co (^,a(/i2; 



Sfifx' 



a{(f_ 



27 
116 



(1 + 2.1930 a{q^) + 6.1460 a^{q^) + 20.5477 c?{q^)) 



X 



(1 + 0.0208 a(/i2) + 0.0095 o?{^l^) + 0.0164 a^(/i^)) 



(A.6) 



where the loop expansion is given in terms of the MS coupling, a. 

Concerning C2/C0, as it obeys the same differential equation eq. flA.4p . the solution differ from 
the one for the gluon propagator only because of the boundary condition which will be now 
obtained from eq. (9) of ref. [20]- Thus, one obtains at the four- loop order: 



C2 l^^,a(/i2; 
Co (^,a(/i2; 



■2/,.2^ 



?,g (/i 



a(g2 



27 
116 



(1 + 1.1728 a{q^) + 2.7098 a^{q^) + 8.4690 oi\q^)) 

(A.7) 



X (1 + 0.0208 a(/i2) + 0.0095 a2(/x2) + 0.0164 a^(/i2)) . 
Thus, we can combine both eqs. flA.6|A.7p . as done in eq. fl2.19p . to obtain: 

. 2^ pert. 2^ A ^ 9 ^VK^^W /^ "('Z') \ ^ 

X (1 + 1.5129 a(g2) + 3.8552 -d^[c^) + 12.4952 a^(g^)) 



X 



(1 + 0.0208 a(/i2) + 0.0095 a^(/i2) + 0.0164 a^(/i2)) . (A.8) 



Finally, provided that the purely perturbative part of the OPE formula in eq. (1A.8P is the coupling 
renormalized in the T-scheme, it seems more appropriate to expand in terms of a^ instead of a. 
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Then, one can apply eqs. fl2.9f2.10p for the conversion and obtain: 

X 



27 
116 



g2 4(iV2-l) V«T(/i2) 

(l + 1.2932 ariq^) + 1.9363 ^^(g^) + 3.8296 a^q^)) 

(1 - 0.7033 ar(/i^) - 0.3652 a^/i^) + 0.0051 ^^.(/i^)) 



.(A.9) 
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